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ABSTRACT 

This paper is aimed at providing a simple model for determination of critical loads, for simple sway frames. It 
is a mathematical model, which is derived systematically from a simplified elastic linear analysis, combined 
with moment equilibrium at the base of frame at the onset of a mechanism. The deflection due to direct axial 
compression and the deflection due to shear deformation are neglected to obtain the overall stiffness. Results 
obtained from the proposed model are striking. 
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I. INTRODUCTION 

The determination of the critical loads for sway frames is associated with complicated mathematics. 
Thus, engineers have tended to develop simplified methods and models for the calculation of critical loads. 
Indeed, the various design manuals are replete with approximate formulae for determination of critical loads. 
The investigation of approximate or simplified models for critical loads has received revived attention in the 
works of Stevens 1 , Home 2 , Anderson 3 , Wood 4 , Goldberg 5 , Orumu 6 , and Hoenderkamp 7 , among several others. 
The basic approach in these solutions is common. It consists in the assumption of a mechanism at the feet of 
the frame (Fig.l). 
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Fig. 1(a) Sway Frame before Critical State 



Fig. 1 (b) Sway Frame at Critical State 



Wheren = ^ j=^ 



The restoring force from sides sway is then determined from the equation of limiting equilibrium about the 

feet. 
QL - PA B = KLA B 
Therefore, 



A„ = 



QL 



fKL-Pl 



At critical stability condition (instability), P reaches its value Per as displacement i B becomes infinitely large. 

In order words, the denominator, in the limit, becomes zero. 
Thus: 
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KL - Per = 
Per = KL 1 



The major difference in the approaches of various researchers consists in the determination of 
the frame structural stiffness k. In this regard, Stevens 3 and Home 8 have applied the rigorous linear 
elastic analysis, which again involves computer application, repeatable for each solution. On the other 
hand, Bolton 9 has determined the critical load for frames working from the tabulated no-shear stability 
functions based on the works of Livesley and Chandler 10 . However, for reasons of their complexity, 
simple formulae that can be applied to any given frame could not be arrived at. 

The author of this paper has derived the stiffness K from a simplified linear elastic analysis 
akin to the approaches by Stevens and Home. However, the proposed simplified method enables the 
stiffness K to be obtained algebraically as a function of the relative stiffness of columns and beams. 

II. THE PROPOSED APPROXIMATE ELASTIC MODEL 

The adopted approach is to determine the sway i B due to a unit horizontal load Q. The stiffness 
K is the reciprocal of i B for the unit horizontal load. For simplification, the effect of shear deformation 
and direct axial compression are neglected. These deflections constitute on the average, about 10% of 
the total deflections 11 . With this assumption, the element slope deflection equation reduces to: 

M AB = k[40 A +20 B ] (2) 

M BA = k[20 A +0 B ] (3) 

The equilibrium of the frame in structure co-ordinates can be written as: 

{P} = {K s } {5} (4) 

Where 

{K s } = structure stiffness matrix; {P} = load vector and {5} = displacement vector. 

For a unit lateral load (Q = 1) applied at joint B (Fig. lb) the proposed simplifies elastic linear analysis to be 
described subsequently yields 

A R = r 

Where 

L = height of frame, Kj = stiffness of column 1 and Y is a dimensionless parameter which depends on relative 
beam to column stiffness ratios in the frame. 

2.1 The Algebraic Approach 

2.1.1 Frame with Rigid Supports 

Consider the frame with relative stiffness n and j as defined in Fig. 1 Let d m and P m represent 
global (master) freedoms and corresponding forces while d s and P s represent local member (slave) 
freedoms and corresponding forces. If member forces and freedoms are related through member stiffness 
K s as follows: 

P s = K s d s 

The member freedoms are related to the master freedoms for compatibility through the transformation matrix 
T as follows 3 

D s =Td m 

Then using the principle of contragradient transformation. 

P m = T T K s Td m 

Thus, the structure stiffness becomes 

K = TtIQT 

The moment equilibrium equation (4) can be solved using the following procedure 

Step 1 : Find moments in elements axes: P s = K s d s 
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Step 2: Determine the global degree of freedom DOF 

Step 3: Set up the system compatibility equation: 
d m = Td m 
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The global stiffness matrix (K = T K S T) after some matrix manipulation, becomes 



K= 



6(I+j) L- -3 L -3j L 

2ki-3 L 2(l-n) n 



V^ 



3jL 



n 2(n+j; 



Where n and j are relative stiffness defined at Fig. lb. 

Step 4: Establish the overall equilibrium matrix for the frame using (1,0, 0) T as load vector. 



= 2k! 



6(H) L- -j L 
- 3 L 2(l+i») 

-3j L n 



' JL 1 
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<*B 
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Step 5: Determine i B by any appropriate method of solution of simultaneous equations. Using method of the 
elimination, A B has the expression. 



1 6Kl 



2(l+j) 



l.5j' L5J Ctl+2;1l[2C■^l+j^j-j^f1;^ 
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Where Y denotes the expression in expression in figure bracket Global frame stiffness equals. 

K=fr 6 



back substituting eqn. 1 the critical load is determined from 



6Kl <„ 
Per = — Y 



Pcr = ^Y 



Where n 2 P E = Cripping load of the column 

Now define a load factor A.crs, such that 

P cr = tars (P B + P c ) 

Per = Pb (1 + a)A.crs 

Where a represents the ratio Pc/Pb 

Then from eqn. 7 

6YPE 

A.crs = — 7— — , 8 

PB i,l+c0ll a 

2.1.2 Frames with Pinned Supports 

If sway frame is pinned at the base, the simplified elastic linear analysis yields the following 
expression for Y 

t = ^4n + 3j + 4nj- 'y-?' , 9 

For symmetrical frame 

If j=l 

Y= D 



-— 8n + 3-7 rr r 7 10 



[ion 
A,crs is determined from equation 8 by substitution of relevant expression for Y 

2.2 Determination of Effective Length 

The effective length of the entire frame can be obtained from the known expression. 

T _ T ' HE. 

Leff — L,\ " 

Per 
Where P cr is the critical load obtained. 

2.3 Modification of model 

2.3.1 Portal frame with fixed base 

The value of Y from the application for frames rigidly fixed at base varies from 1 to 4 for 
KB/KC ratio of to oo respectively. The expected value of A.cr for this range is 0.25 to 1. This is 

f Y= 

equivalent to the range -^— to — — which means the critical load factor earlier obtained for every 

KB/KC needs to be multiplied by the ratio of 0.25 and further multiplied by the amplification factor f amp 
which was derived based on the maximum deflection for KB/KC=0 which is the case of a cantilever. 

f amp = 1 + 0. 083333 Sin ^^^ 11 

The modified model becomes 

A.crs = ^ PE , xl.214876 fa™ 

P5 ,'--z : -!- p 

D.739Y PE 
Xcrs = ' f^p 12 

PS (1+1} 
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2.3.2 Portal frame with pinned base 

The value of Y from the application for frames pinned at base varies from to 1 for KB/KC 
ratio of to oo respectively. The expected value of Xcr for this range is to 0.25, which means the 
critical load factor earlier obtained for every KB/KC needs to be multiplied by the ratio of 0.25 and 
further multiplied by the amplification factor f amp which was derived based on the maximum deflection 
for KB/KC=0 which is the case of a cantilever. 

f amp = 1 + 0. 083333 SirmY 13 

The modified model becomes 



tars = ^ PE , xl.214876 f a , 

P5 '.-Z.-I- 



. D.7J9\ PE. 

A.crs = — fai 

PE [H-a] 



14 



III. APPLICATION 1 

Find the critical load factor and effective length for various values of n (K B /K C ) for the frame 
shown in fig. la, where j = 1, if the frame is fixed at base and if the frame is pinned at base. 
Solution: 
Substituting various value of n for j=l, the results are shown in table 3 and table 4 below will be obtained 

Table 3: Values of elastic critical load factor (A,cr) for Sway Frame with fixed feet. 



n 


f 


Acr 


Acr 
irodifed 


Acr 
Home 


Kdiff 


effective 
length 





1 


0.303719 


0.25 


0.25 





2 


0.1 


1.391 


0.422565 


0.35937394 


0.36 


-0.174 


1.668117 


0.25 


1.B1B 


0.552216 


0.48317233 






1.43863 


0.5 


2.286 


0.694215 


0.61785371 


0.61 


0.129 


1.272205 


0.75 


2.588 


0.786096 


0.70075037 






1.194589 


0.8 


2.535 


Q.S00714 


Q.7134551 


0.7 


1.922 


1. 183904 


0.85 


2.681 


0.S14367 


0.72518652 






1.174289 


0.9 


2.723 


0.S27149 


0.73604304 






1.165596 


0.95 


2.763 


0.839141 


0.74611437 






1.157703 


1 


2.8 


0.850413 


0.7554783 


0.75 


0.730 


1.150506 


1.5 


3.077 


0.934520 


0.82198614 






1.102980 


2 


3.25 


0.9S70S7 


0.86037702 






1.078091 


4 


3.571 


1.0S4711 


0.92514016 


0.94 


-1.581 


1.039672 


5 


3.647 


1.107681 


0.93921199 






1.031854 


8 


3.769 


1.144787 


0.96110011 


0.97 


-0.918 


1.020036 


9 


3.793 


1.152038 


0.96526347 






1.017834 


10 


3.813 


1.157929 


0.96862046 


0.98 


-1.161 


1.016069 


BO 


4 


1214815 


1 


1 





1 
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Table 4: Values of elastic critical load factor (A.cr) for Sway Frame with pinned feet 



n 


T 


Acr 


Acr 
nrodified 


Acr 
Home 


% 
diff 


effective 
length 




















= 


0.1 


0.57 


0.0505 


0.043402 


0.04 


8.5 


4.80000 


0.2 


0.285 


0.0B67 


0.075082 


0.075 


1.4 


3.62541 


0.25 


0.5 


0.1518 


0.13541 


O 




2.71746 


0.5 


0.6 


0.1822 


0.161888 


O 




2.48537 


0.75 


0.615 


0.1869 


0.165833 







2.45563 


0.E 


0.63 


0.1912 


0.169452 


0.167 


1.46E 


2.42927 


0.85 


0.643 


0.1952 


0.172780 







2.40575 


0.9 


0.666 


0.2024 


0.178695 







2.36561 


1 


0.75 


0.2277 


0.198548 


0.184 


7.907 


2.24423 


1.5 


0.8 


0.2429 


0.209795 


O 




2.18324 


2 


0.889 


0.2699 


0.228556 







2.09172 


4 


0.909 


0.2761 


0.232608 


0.233 


0.168 


2.07342 


5 


0.941 


0.2858 


0.238897 


O 




2.04594 


8 


0.947 


0.2877 


0.240090 


0.242 


0.788 


2.O4085 


9 


0.952 


0.2892 


0.241052 







2.03678 


10 


0.999 


0.3035 


0.249907 




2.421 


2.00036 


« 


1 


0.3037 


0.25 


0.25 





2 



3.1 Application 2 

Find the critical load factor and effective length for various values of n (K B /K C ) for the frame 
shown in fig. 2, for unsymmetrical portal frame i.e j = to oo, if the frame is fixed at base. 



Solution: 

Substituting various value of n for j=n, the results are shown in table 5 and table 6 below will be obtained 



XI 



- 1 - Fig2i 
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Table 5: Values of elastic critical load factor (A.cr) for Sway Frame with fixed feet. Unsymetrical j=n 



n&j 


f 


Acr 


kr 

irodifed 


effective 
length 


Acr 

Home 


Kdiff 





2 


0.6074 


0.536084 


1.36578 






0.1 


2.0241 


0.6147 


0.543075 


1.35696 






0.2 


2.0580 


0.6281 


0.555774 


1.34137 






0.25 


2.0961 


0.6366 


0.563861 


1.33172 






0.5 


2.2826 


0.6932 


0.51698 


1.27310 






0.75 


2.5235 


0.7664 


0.683455 


1.20960 






1 


2.8 


0.8504 


0.755478 


1.15050 


0.75 


-0.730 


2 


4.0909 


1.2424 


1.01462 


0.99276 






4 


7 


2.1260 


1.75 


0.75592 






5 


8.51 


2.58 


2.305219 


0.55863 






10 


16.2 


4.93 


3.97677 


0.50145 


4 


-0.580 


10000 


1571. 


4773 


3805.86 


0.01621 
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Table 6: Values of elastic critical load factor (tar) for Sway Frame with fixed feet. Unsymetrical j^n 



n 


i 


T 


J.ir 


Ae r m □ d if s d 


Alt H = -i 


= ff^£tive length 


■3 


1 


1 


■D.3D37 


■D.2E 


■3.2 5 


2 


■D.l 


1 


1.3913 


0.4225 


■3.3 59 37 


•3.3 6 


1.66B1 


■D.l 


2.5 


■D . 2D 3 6 


O.D633 


■D.D4 3 96 




4. 5193 


■3.1 


5 


-1.776 


-3.539 


-3.4 3 54 5 




- 


■D.l 


10 


-E.7S2 


-1.747 


-1.35332 




- 


■3.25 


1 


1.3132 


-3.5522 


■3.4 3 317 




1.4336 


0.25 


2.5 


1.3411 


0.3162 


■3. 2612-37 




1.95663 


•3.2 5 


5 


-3.315 


-3.095 


-3.37232 




- 


■3.2E 


10 


-3.05S 


-3 . 9 2.3 


-3. 32-359 




- 


■3.5 


1 


2.2357 1 


•3.6942 


•3.617B54 


■3.61 


1.27221 


■3.5 


2.5 


2.0352 


•3 . 6 13 1 


■3.546279 




1.3529 3 


■3.5 


5 


1.4656 


0.4451 


3 . 3 S3 7 1 B 




1.62-363 


■3.5 


10 


•3.2470 


0.0750 


■3 .-3 5 B 1-3 3 




4.1435 9 


1 


1 


2.B 


•3. 3 5-34 


•3.75547B 


O. 75 


1.15-351 


1 


2.5 


3.25 


0.9370 


'3.36-3377 




l.-373'3 9 


1 


5 


3.7021 


1.1244 


•3.94 9 2-3 2 




1.-32641 


1 


10 


4.413H 


1. 34-36 


1.-364B3B 




-3.969'3S 


2 


1 


3.25 


.937 1 


■3. B 6-3 377 




1.-37B-3 9 


2 


2.5 


4.45 


1.352 


1. -3 7-34 11 




'3.96655 


2 


5 


■6 


1.3223 


1.391747 




•D . 34 7 6 5 


2 


10 


H.7143 


2.6467 


2.3ESE67 




•D. 65155 


5 


1 


3.6471 


1.1077 


'3.939212 


■3.94 


1 .D 3 1 3 5 


5 


2.5 


5.6634 


1.7219 


1.3-3 I'D S 




'D. 37669 


5 


5 


3.5113 


2.53E1 


2.30E219 




■D.6ES63 


5 


10 


13.4 3 3 


4.0979 


3.E1177E 




•D . E 3 3 6 2 


10 


1 


3.3125 


1.1579 


■3.96362 


•3.9 3 


l.'D16-37 


10 


2.5 


6.25 


1.3932 


1.47-34 2 9 




■D . 3 24 6 7 


10 


5 


9.3571 


2.9935 


2.494393 




■3. 63310 


10 


10 


16.231 


4.929S 


3.97677 4.3 


7.E17 



IV. DISCUSSION 

The results from the worked examples demonstrate the simplicity and accuracy of the proposed 
mathematical models for simple sway frames. The solutions from Home and Merchant (2) using the 
rigorous matrix computer analysis using stability functions are tabulated and used as accurate baseline 
values for comparison with the critical load and modified critical load factors. The effective lengths for 
each case are also presented. Tables 3,4,5 and 6 show that the results are very reliable when compared. 
A rather controversial problem is discussed below For a particular case where the beam is 101, a stocky 
column is 101 and a slender column of the sway frame is I, four buckling loads are predicted by the 
model viz. 
From table 6 for 

n=l, j=10 tar modlM = 1.065 
n=10, j=0.10 Icr mo difi e d=0.59 
n=l, j=0.101cr modlfled =0.57 
n=10,j=10 tar modlM =3.98 



101 



1 01 



The two of these critical load factors to be given attention to are the once with integers and not fraction for j i.e 
n=l, j=10 tar modlfied = 1.065 
n=10,j=10 tar modlfied =3.98 

Wood 10 had obtained a value of tar =4.3 and several discussants of his paper had argued that this value 
cannot be possibly the critical load and claimed that the computer program may have skipped the 
fundamental critical load which according to the discussants should be around the critical load factor of the 
slender column which for the KB/KC ratio of 10 is in the region of tar =0.98 or slightly above because of the 
support it will receive from the stocky column. The results obtained in the method herein discussed have 
validated their arguments. The critical loads Per =1.065PE and Per =3.98PE appear to be the first and second 
mode buckling loads of the special frame. The frame will have failed before hitting Per =3.98PE speaking 
from Engineering point of view. Therefore the critical load is Per =1.065PE. 

In deriving the analytical models, deformations due to axial and shear stress resultants have been 
ignored. As practice has shown, these components contribute for less than 10 percent of the total stresses in 
the systems stability response. This, perhaps, also explains the minor disparity between the results from the 
proposed models and those obtained from the rigorous analytical efforts of Home and Merchant. 
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CONCLUSION 

The proposed analytical models offer a simplified hand of solution to sway stability problems of 
building frames. Its success is justified by close reproduction of existing solutions of rigorous methods. Its 
application will have tremendous impact on reducing cost of computing and increasing scope of solvable 
problems on sway stability of frames, while opening a greater scope for studies in this rather complex area of 
structural analysis and design. 
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